Terahertz surface plasmon polariton propagation and focusing on periodically 

corrugated metal wires 
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In this letter we show how the dispersion relation of surface plasmon polaritons (SPPs) propa- 
gating along a perfectly conducting wire can be tailored by corrugating its surface with a periodic 
array of radial grooves. In this way, highly localized SPPs can be sustained in the terahertz region 
of the electromagnetic spectrum. Importantly, the propagation characteristics of these spoof SPPs 
can be controlled by the surface geometry, opening the way to important applications such as energy 
concentration on cylindrical wires and superfocusing using conical structures. 
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One of the major driving forces in the field of plas- 
monics is the ability to spatially confine electromagnetic 
energy at visible frequencies over distances significantly 
smaller than the wavelength. Surface plasmon polari- 
tons (SPPs) and localized plasmons at metal/dielectric 
interfaces open up a previously inaccessible length scale 
for optical research, with promising applications in the 
miniaturization of photonic circuits Q, Q , near-field op- 
tics Q and single- molecule optical sensing 0,0. It would 
be greatly advantageous to take concepts such as hi ghly 
localized waveguiding [j| 0, and superfocusing IfTllfl 
to lower frequencies, particularly the THz regime 1 1 ill - 
where plasmonics could enable near-field imaging |l2( 
and biosensing 13] with unprecedented sensitivity. How- 
ever, there is a basic problem: SPPs offer sub- wavelength 
field localization only for frequencies close to the intrin- 
sic plasma frequency of the conductor, which for most 
metals is in the ultraviolet part of the spectrum. At 
THz frequencies on the other hand, metals resemble in 
many ways a perfect conductor, and the negligible pene- 
tration of the electromagnetic fields leads to highly delo- 
calized SPPs akin to grazing-incidence light fields. In this 
frequency range, SPPs are also known as Sommcrfeld- 
Zenneck waves |l4j . Such weakly guided THz waves have 
recently been characterized on cylindrical metallic wires 
(la IWL Il7| . While this waveguide geometry is in princi- 
ple suitable for a variety of applications such as biologi- 
cal sensin g^ in confined spaces jl5j or the local probing of 
materials [lg, the radial extent of the fields over many 
wavelengths and the concomitant radiation loss due to 
bends, non-uniformities or nearby objects severely limit 
the efficiency and viability of such concepts. 

In this letter we present a very promising route to 
achieving sub-wavelength confinement in metal wires at 
THz or microwave frequencies by showing that the dis- 
persion relation of SPPs can be engineered at will by peri- 
odically structuring the cylindrical surface with grooves. 



In this way, the confinement does not rely on the finite 
conductivity of the wire, but is purely due to the sur- 
face structure. It is interesting to note that the viabil- 
ity of this approach was hinted at by Goubau in 1950 
|14| . More recently, the idea of tailoring the topography 
of a perfect conductor to allow the existence of surface 
waves resembling the behavior of SPPs at optical fre- 
quencies was discussed in the context of two-dimensional 
hole lattices 0] and one-dimensional groove arrays [2(J 
machined into flat interfaces. Experimental verification 
of these findings in the microwave regime has been re- 
cently reported 0| . Due to their mimicking characteris- 
tics, these geometry-controlled surface waves were named 
spoof SPPs. Here we demonstrate that spoof SPPs can be 
sustained and exhibit a rich behavior in three dimensions, 
including the possibility of deep sub-wavelength energy 
concentration in cylindrical and superfocusing in conical 
structures. The length scales of the surface structures 
that we describe are within the limits of state-of-the-art 
laser-machining techniques that are used to make corru- 
gated feedhorns for sub-millimeter radiation [i^ . 

The inset of Fig. la shows a schematic picture of the 
envisioned structure. An array of radial grooves of depth 
h = R—r, width a and lattice constant d is machined into 
a perfectly conducting cylinder of radius R. Within the 
perfect conductor approximation, valid for metals in the 
microwave or THz ranges of the electromagnetic spec- 
trum, the frequencies of the dispersion bands scale with 
the reciprocal size of the structure. Therefore, in Fig. 1 
d is used as the unit length. We are interested in calcu- 
lating the dispersion relation (frequency, u>, versus axial 
propagation constant, k) of TM-polarized waves propa- 
gating in the z-direction along the wire. First we apply 
a theoretical framework based on a modal expansion of 
the EM-fields. In order to yield quasi-analytical expres- 
sions for io(k), we search for surface waves that have no 
azimuthal (<fi) variation, i.e. the azimuthal order m is 



equal to zero, and when describing the EM-fields inside 
the subwavelength grooves, only the lowest eigenmode is 
considered in the expansion. In this way, in region II, 
E and B are zero everywhere except inside the radial 
grooves and E z can be written as: 



E I z I (p)=AJ (gp) + BN (gp) 



(1) 



where A and B are constants, Jo and Nq are the zero- 
order Bessel and Neumann functions, respectively, and 
g = lu/c, with c being the velocity of light. Within this 
single mode approximation E 1 ^ = 0. 

On the other hand, the z-component of the E-field in 
region I (vacuum) can be expressed as: 



+00 



(2) 



where C n are constants, k n = k + 2nn/d takes into ac- 
count diffraction effects and q n = yfe^ — g 2 . The radial 
dependence is controlled by Kq, the zero-order modified 
Neumann function, which displays a non-divergent be- 
havior as p — + 00. Note that, as we are interested in the 
non-radiative region of the dispersion relation, k must 
be outside the light cone, i.e., k > g. The other non- 
zero components of E and B can be obtained straight- 
forwardly from E z . The dispersion relation of the TM- 
polarized waves propagating along the cylinder is deter- 
mined by imposing matching conditions on the parallel 
components of E and B at the interface between regions 
I and II and at the bottom of the grooves. Following this 
procedure yields a transcendental equation for u)(k): 



+00 

E 



2 g K^R) N Q (gR)J Q (gr) - N (gr)J (gR) 
n q n K Q {q n R) N (gr)Ji(gR) - N 1 (gR)J (gr) 



= 1 



(3) 

where Ji and N± are the first-order Bessel and Neu- 
mann functions, respectively, and K\ is the first-order 
modified Neumann function. S n is given by S n — 
y/ajd sinc(fc n a/2). 

Figure la shows the dispersion relation calculated us- 
ing Eq.(3) for SPPs propagating along a periodically cor- 
rugated, perfectly conducting cylinder. The inclusion of 
just 3 or 4 diffraction orders in the sum appearing in 
Eq.(3) yields sufficient convergence. The width of the 
grooves is a = 0.2c? and two different depths of inden- 
tation are analyzed: h = R — r = QAd (dashed lines) 
and h = d (full lines). For both depths, several values 
of the outer (R) and inner (r) radii are shown for fixed 
h. Figure la clearly shows that the inclusion of a peri- 
odic array of grooves on a perfectly conducting wire dra- 
matically changes the dispersion relation of SPPs. Note 
that, within the perfect conductor approximation, SPP 
bands coincide with the light line for the non-corrugated 
geometry, and the modes are not bound. For the case 
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FIG. 1: (a) uj versus k for the corrugated, azimuthally invari- 
ant wire geometry shown in the inset for two different cases: 
h — R — r = 0.4d (dashed lines) and h = d (full lines). For 
each value of h, the dispersion relations calculated for increas- 
ing values of R and r (h fixed) are presented. The dotted line 
shows the band obtained with a FDTD code for h = d and 
R = 2d. (b) Dispersion relation of the SPP bands obtained 
with the FDTD code for azimuthal modes m = 0, m = 2 and 
m = 4. The insets show the amplitude of the E-field in a 
quarter cylinder for these three bands at k — n/d. 



of spoof SPPs, ui{k) departs significantly from the light 
line and this departure is greater when the indentations 
are deeper. This dispersion relation closely resembles 
the behavior of SPPs propagating along a metal wire at 
optical frequencies. In this last case, for m — 0, u(k) ap- 
proaches uo p /\/2 (where u> p is the plasma frequency of the 
metal) at large k. For the spoof SPPs bands, the asymp- 
totic frequency (u>s) is mainly controlled by the depth 
of the grooves. An analytical expression for ws can be 
obtained in the limit R,r » d and A ^> d, a. By us- 
ing the asymptotic expansions of the different Bessel and 
Neumann functions appearing in Eq.(3) and neglecting 
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the diffraction orders, k(u>) becomes: 



k = gJl + -pt&n 2 gh (4) 

The asymptotic frequency ujs — nc/2h is inversely pro- 
portional to h, explaining the behavior of u(k) observed 
in Fig. la when the depth of the radial grooves is varied. 

In order to check the accuracy of the approximations 
made in our quasi-analytical model, we have conducted 
finite-difference-time-domain (FDTD) numerical simula- 
tions on one unit cell of the wire structure. In Fig. la 
the spoof SPP band obtained with this code for the case 
of h = d and R = 2d is also plotted (dotted line). The 
agreement between the band obtained from Eq.(3) and 
the FDTD result is remarkable. In Fig. lb we present 
the to — 0,2,4 bands (only for k > g) for the corru- 
gated wire analyzed above calculated with the FDTD 
code. The spatial dependence of the amplitude of the 
E-field for the three bands (to = 0,2,4) evaluated at 
k = 7r/d, border of the first Brillouin Zone, is plotted 
in the inset. The highly localized character of the spoof 
SPPs along the radial direction is evident. 

Once it has been shown how the bands of spoof SPPs 
can be engineered, we now focus on their light guiding 
properties in a wire tailored to operate at a frequency 
near 1 THz. In panel (a) of Fig. 2, we show the disper- 
sion relation of the spoof SPPs bands obtained with our 
quasi-analytical approach for a — 10 /im, d — 50 /j,m, 
r = 50 fim and R = 100 fim (black line). For this 
set of geometrical parameters the asymptotic frequency, 
fs, is equal to 1.1 THz. The spatial variation of the 
E-field associated to the spoof SPPs for different fre- 
quencies is displayed in panels (b-d) of Fig. 2. These 
pictures have been obtained by numerical simulations 
using the finite integration technique (FIT). The wire 
has been excited using an input port with a radially po- 
larized broadband THz pulse, and multiple calculations 
with adaptive mesh refinement ensured well-converged 
results. For / = 0.6 THz, the E-field is highly concen- 
trated on the surface of the wire but this concentration 
is much stronger at a higher frequency, / = 1.0 THz. 
This is due to the fact that in this last case the prop- 
agation constant is much longer than u/c (see Fig. 2a) 
resulting in a very large radial component of the wavevec- 
tor, q = yjk 2 - g 2 . For / = 1.2 THz (> f s ), no spoof 
SPPs propagate along the corrugated wire and light is 
scattered at the input port. 

Both the high localization and low group velocity of 
spoof SPPs on corrugated metal wires are advantageous 
for applications in the routing of radiation with low bend 
losses and for sensing. The fact that the dispersion rela- 
tion of the spoof SPPs can be easily tuned by engineering 
the geometry of the indentations opens up an even more 
fascinating possibility, the gradual energy concentration 




FIG. 2: (a) Dispersion relations of spoof SPPs for a = 10/im, 
d = 50/xm and R — 100/xm for five different values of the inner 
radii r = 85, 60, 50, 40 and 20/im. (b-d) Spatial variation 
of the E-field amplitude (logarithmic scale) calculated using 
the finite integration technique on a wire with a = 10/im, 
d — 50fj,m, R — 100/xm, and r = 50[im for three different 
frequencies, / = 0.6, 1.0 and 1.2 THz. 



of radiation propagating along wires, and superfocusing 
of radiation at the tip of conical structures. 

Panel (a) of Fig. 2 gives us a clue on how THz radi- 
ation can be concentrated using periodically corrugated 
metal wires. Let us imagine a metal wire in which the 
outer radius is fixed (R = 100 jim) and the inner one 
(r) is adiabatically (i.e., back-reflection and scattering 
are negligible) reduced, thus increasing the groove depth, 
from r = 100 fim at the beginning to 20 fim at the end of 
the wire. If the frequency of the surface wave is fixed (for 
example, / = 0.6 THz, A = 500 fim, see dotted line in 
Fig. 2a), it is expected that the group velocity of the mode 
will be gradually reduced as it propagates along the wire 
in the direction of decreasing r. Along with this change 
in the group velocity, the spatial extension of the mode in 
the radial direction would be severely reduced, as can be 
seen in Fig. 3 (dotted lines), in which the dependence of 
the E-field concentration with r is analyzed. Notice that 
these curves have been obtained with the quasi-analytical 
method, assuming that the wire is infinite and r is fixed. 
To confirm the idea of gradual energy concentration, we 
have conducted FIT numerical simulations on a finite 
straight wire of radius R = 100 fim and length 2 mm 
containing 40 grooves with d = 50 fim. In this wire the 
inner radius r decreases from 95 fim to 20 fim in steps of 
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FIG. 4: Superfocusing on a corrugated cone of length 2 mm, 
with constant h — 5 /im and d = 50/im. 7? is reduced from 
100/ira to Plot shows the magnitude of the E-field on 

a logarithmic scale spanning two orders of magnitude. 



FIG. 3: Field concentration via adiabatic reduction of r. Dis- 
tribution of the E-field, evaluated at / = 0.6 THz, along the 
radial direction at different locations of the wire correspond- 
ing to r = 85, 60, 40 and 20/j,m calculated using the quasi- 
analytical method for an infinite wire (dotted lines) and FIT 
for the real finite structure (solid lines) on a linear color scale. 
The inset shows the distribution of the E-field in a logarithmic 
scale. 



5 /im every two lattice periods. This region is followed 
by a number of lattice periods with constant r — 20 fim. 
The distribution of the E-field for / = 0.6 THz is plotted 
on a logarithmic scale in the inset of Fig. 3. As expected 
from our previous discussion, the mode becomes increas- 
ingly confined to the surface with propagation distance, 
as shown in the main part of Fig. 3 (solid lines). At the 
end of the wire, the E-field is concentrated on a region 
10 times smaller than the wavelength on the surface of 
the wire. The agreement between the quasi-analytical 
method and FIT validates the idea of concentrating light 
via an adiabatic reduction of r. 

The same concept can be extended to conical struc- 
tures, where both R and r are adiabatically reduced. In 
this case, we expect focusing as long as the gradual re- 
duction in r is greater or equal to that in R. Figure 
4 demonstrates this effect for a cone where the groove 
depth R— r = 5 [im is constant, d = 50 [im and R is re- 
duced from 100 [im to 10 [im over a distance of 2 mm. As 
apparent from this figure, efficient deep sub-wavelength 
confinement to the tip (2r = 0.04A) takes place. This ge- 
ometry has great promise for channeling THz radiation to 
micron-scale volumes for near-field imaging, spectroscopy 
and sensing applications. 

In conclusion, we have shown that periodically cor- 
rugated metal wires can sustain spoof surface plasmon 
polaritons, and that their dispersion and mode profile is 
determined by the geometry. Apart from low-loss propa- 
gation due to the good confinement, such structures allow 



energy concentration and superfocusing. We expect these 
results to enable a new class of THz optical research and 
technology with only minor extensions of existing micro- 
fabrication techniques. 
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